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We propose a three terminal heat engine based on semiconductor superlattices for energy har¬ 
vesting. The periodicity of the super lattice structure creates an energy miniband, giving an energy 
window for allowed electron transport. We find that this device delivers a large power, nearly twice 
than the heat engine based on quantum wells, with a small reduction of efficiency. This engine also 
works as a refrigerator in a different regime of the system’s parameters. The thermoelectric perfor¬ 
mance of the refrigerator is analyzed, including the cooling power and coefficient of performance in 
the optimized condition. We also calculate phonon heat current through the system, and explore the 
reduction of phonon heat current compared to the bulk material. The direct phonon heat current 
is negligible at low temperatures, but dominates over the electronic at room temperature and we 
discuss ways to reduce it. 


I. INTRODUCTION 

There has been increasing interest in developing high 
efficiency, high power thermoelectric devices, constructed 
from the bottom-up using nanoscale designs. The pri¬ 
mary applications driving interest in this area are energy¬ 
harvesting, the collection and conversion of waste heat 
to electrical power, produced from sources ranging from 
hand-held electronics to industrial sources of heat, and 
refrigeration, actively cooling a spatial region via elec¬ 
trons to evacuate heat out of an area. The use of 
nanoscale architecture instead of bulk materials is moti¬ 
vated by the low figure of merit - or poor thermoelectric 
conversion efficiency - of bulk materials, whereas con¬ 
duction through nanoscale electronics can reach Carnot 
efficiency. 

One way to produce high thermodynamic efficiency in 
the conversion of heat to power is the use of structures 
with sharp spectral features, such as quantum dots [D- 
[3] . The use of quantum dots in thermoelectric transport 
has been extensively researched in the past several years 
[4HTT] . See Ref. [1^ |T3] for recent reviews of these and 
related activities. In particular, if the dot is transporting 
electrons via resonant tunneling, the quantum dot acts 
as an energy filter - permitting the “tight-coupling” of 
heat and charge transport which can lead to Carnot effi¬ 
ciency. Other structures from mesoscopic physics, includ¬ 
ing the quantum point contact and electron cavity [14], 
quantum wells [m HE], quantum Hall bar [m ng , su¬ 
perconducting leads m, and Coulomb blockaded quan¬ 
tum dots(s) [2Qll24] have also been investigated for their 
multi-terminal thermoelectric properties. The late Prof. 
Markus Biittiker, for whom this special issue is in mem¬ 
ory of, was highly influential in the theoretical develop¬ 
ment of these ideas, as can be seen in the above list of 
references. 
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Several experiments have begun exploring this physics. 
Prance et al. [25] performed experiments on a cavity con¬ 
nected to resonant tunneling quantum dots acting as an 
electronic refrigerator, based on the proposal of Edwards 
et al. m- They demonstrated that applying bias to 
the system results in cooling a large Q/irin? cavity from 
280mi^ to below IQOmiC. Very recently, Roche et al 
m and F. Hartmann et al. [28] showed rectification of 
electrical current of the nano Amp scale and power pro¬ 
duction on a pico Watt scale from a capacitively coupled 
source of fluctuations. This was based on the theoretical 
proposal of Sothmann et al. M- 

While a nanoscale thermoelectric generator can power 
nanoscale devices, it is of great interest to find practical 
ways to scale up these nanoengines. One way is to simply 
add them in electrical series while being able to couple to 
a common source of heat. In commercial thermoelectric 
generators, this is usually done by alternating the semi¬ 
conductor type, of either p-type or n-type to be able to 
apply the heat difference in parallel because the heat and 
electrical transport are in opposite directions in a p-type 
semiconductor [29]. This permits the generated voltage 
to grow with the number of elements, while keeping the 
current fixed. Various other ways of scaling the devices 
have been proposed [3QH33] . In Jordan et al, a layered 
structure was proposed by alternating layers of semicon¬ 
ductor and self-assembled quantum dots, so as to create 
a large-scale device where heat and electrical transport 
are separated, while keeping the high thermodynamic ef¬ 
ficiency [9]. This is a parallel strategy of scaling, so the 
generated current grows with the number of dots, while 
the voltage difference is fixed. Sothmann et al. consid¬ 
ered a technically simpler method of creating quantum 
wells that permit resonant tunneling m- The physics 
there is somewhat different because energy may be dis¬ 
tributed into the transverse degrees of the electron mo¬ 
tion. Nevertheless, reasonable thermodynamic efficiency 
was found, with increased power production. 

One of the outstanding challenges to creating high- 
efficiency thermoelectric devices is phonon transport. 
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Phonons give a way for the hot and cold side of the de¬ 
vice to exchange energy directly, without converting it 
to power via the electrons. Therefore, any possible way 
to reduce the phonon transport while still allowing elec¬ 
tron transport will aid in the overall thermodynamic effi¬ 
ciency. Interface-based devices, such as described above 
can help with this, because the interface helps to reflect 
the phonons j34H4T| . Ideally, there will be additional ma¬ 
terial layers that act as thermal insulators. 


The purpose of the present article is to build on 
these accomplishments, and make an analysis of a 
thermoelectric device based on semiconductor super¬ 
lattices. These structures are fabricated by making a 
periodic layered structure of alternating materials, such 
as GaAs/AlGaAs. The effect on the electronic transport 
is to form a series of mini-bands of allowed and forbid¬ 
den energies where conduction electrons can transport 
[42II45] . This structure can be considered as a general¬ 
ization of the resonant tunneling quantum wells. The 
mini-band gives a top-hat profile of variable width for 
the energy-filtering. Such a top-hat profile has been ar¬ 
gued by Whitney to offer the highest efficiency for a given 
power extraction [46]. However, we note the transverse 
degrees of freedom make the system somewhat differ¬ 
ent. At a small band width, our system will be sim¬ 
ilar to a quantum well, but can be extended to allow 
a fixed width longitudinal energy window. We make a 
first-principles analysis of the heat and charge transport 
in a three-terminal geometry, where two terminals carry 
charge, and a third carries heat (see Fig. [^. The off¬ 
set of the miniband centers and their respective widths 
determine the power produced and efficiency of the heat 
conversion given fixed temperature differences. We next 
make an analysis of the coefficient of performance of this 
device for the purposes of refrigeration of the central re¬ 
gion. The final purpose of the present work is to also 
make a systematic calculation of the heat transport due 
to phonons. We make a detailed investigate the heat 
current through the device from phonon transport using 
a Kronig-Penney model, and consider different ways to 
stop it. 


The paper is organized as follows. In section [ll] we 
introduce a model of the superlattice heat engine in its 
dual roles: the energy harvester in section II A| and the 
refrigerator in section [XTCj We discuss our results for the 
generated power and the efficiency of the engine in section 
El and show the cooling power and the coefficient of 
performance in section |ll G[ The second part, section O 


focuses on phonon heat current generated by the heat 
engine and discusses the effect on the efficiency of the 
heat engine. We finish with our conclusions in section 

llYl 


II. HEAT ENGINE BASED ON 
SUPERLATTICES 


A. Energy harvesting by electron transport 


We consider a setup shown schematically in Fig. It 
consists of a center cavity connected to two electronic 
reservoirs r = L^R via a superlattice. The electronic 
reservoirs are characterized by the occupation of the 
states given by the Fermi function, fr{E) = [exp[(F^ — 
/ir )/with temperature Tr and chemical po¬ 
tential and the underlying assumption is that inelastic 
scattering processes restore the local thermal equilibrium 
on a fast time scale. The center cavity is also assumed 
to be in thermal equilibrium with a heat bath of tem¬ 
perature Tc. We assume the structure is translationally 
symmetric within the x and y directions, perpendicular 
to the growth direction 2 :. The superlattices are designed 
as periodic structures with lattice constant d (sum of the 
width of well and barrier thickness). The periodicity of 
the structure in the 2 : direction implies that the eigen¬ 
states of the Hamiltonian can be written as Bloch states 
with the Bloch vector q G [—7r/d,7r/d], so the simple so¬ 
lutions with the Kronig-Penney model resemble the stan¬ 
dard superlattices m- The corresponding eigenvalues of 
the Hamiltonian form a miniband. The allowed energies 
of the miniband can be written. 


E^{q) = — 2/3^ cos gd, (1) 

where this is the result of a standard tight-binding calcu¬ 
lation with band indices u. £q is a center of miniband u 
and its width is 4/3^. Our discussion is restricted to the 
electron transport through the lowest miniband, thus the 
miniband index u is neglected from here on. 

To find the electric and heat currents through the su¬ 
perlattices expressed in terms of an integral over energy, 
we first find the density of states of the superlattice. For 
the given miniband, the energy is that of a two dimen¬ 
sional electron gas with the bottom of the band at E{q) 
in Eq. 0 . Therefore the three dimensional density of 
states is given by [47] 

1 

V3d{E) = - de vsL{e) V2 d{E - e) 

^ J — 00 

= ^y deusiie)- (2) 


The factor 1/2 is to avoid double counting the spin and 
U 2 d = vnj'KT? is a two dimensional density of states per 
unit area. The density of states of the superlattice can 
be factorized into longitudinal and transverse parts, and 
the one dimensional superlattice density of states vsl is 


J^sl{E) 


1 Q{E - E:)@{-E + Et) 


( 3 ) 


where we use a Heaviside step function 0 to show the 
range of the energy E with the maximum/miminum value 
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FIG. 1. (Top) Schematic of the superlattice heat engine. A 
hot cavity at temperature Tc is coupled via superlattices to 
cold electronic reservoirs at temperature T^. (Bottom) The 
periodic structure of the superlattices form the miniband cen¬ 
tered at with the width 4/3^^^ when we apply the bias 

voltage [iR — 11L — gV . The gray shading area shows the 
energy miniband where the electrons can transport between 
regions. The shadings in the source, cavity, and drain regions 
indicate thermal smearing. 


Ef = Therefore, the electrons only in selected 

values of energy E~ < E < E^ will transport and gen¬ 
erate current. We assume the central cavity region is 
strongly coupled to the external source of energy, so the 
occupation is described as a Fermi function with local 
temperature Tc determined by the thermal reservoir. 

The electric and energy currents for simplified mini¬ 
band transport in 2 ; direction emitted by the reservoir 
r into an cavity c can be evaluated within a Landauer- 
Biittiker approach as 

dE^dE^lME) - UE)], (4) 

Jr = ^1^ [ dE,dE^ E [fr{E) - f,iE)], (5) 

where A is the surface area of the superlattice, E'^^ is the 
maximum/minimum energy of the reservoir r, and we de¬ 
note E±_ as the energy carried in the transverse degrees 
of freedom, and E^ as the energy carried in the longitu¬ 
dinal degree of freedom, so that E = E± + Ez. Here, the 
square root in Eq. (§ cancels the velocity of the electron, 
Vg, in the current Ir = (e/4) / dE^Vgiy 3 D[fr{E) - fdE)] 
|48j . We see that the range of the integral comes from 
the density of states which gives a transmission function 
of flat box form, T^{Ez) = Q{Ez — El~)&{—Ez + 

To rewrite the above equations, we introduce the inte¬ 
grals Ki{x) = dt{l + _ iog(l + e^) and 

K 2 {x) = = —Li 2 (—e^) with the dilog¬ 
arithm Ij\ 2 {z) = I 2 • Then the simplified analytic 


expressions of the electric and energy currents are 




eu2DAkB 
27rh 

l^2DAkB 


2iTh 


i: 

L 


dE, 


dE,E 


TrKdm - TMEd\ , (6) 

TrKdm - TcKdm 


/' 

JEl 


-\-kB I dEz 


T^KdEl) - 


, ( 7 ) 


where E^ = (/i^ — Ez)/kBTr and E^ = (/ic — Ez)/kBTc 


B. Results 


We now analyze the system by focusing on linear re¬ 
sponse and later turn to the nonlinear regime. To sim¬ 
plify the analysis of the system, we introduce the aver¬ 
age temperature T = (T^ + ^c)/2 and the temperature 
difference AT = Tc — T^. For energy harvesting, the 
temperature difference is considered to be AT > 0. (If 
we consider refrigeration, the temperature difference is 
AT < 0.) We introduce the bias /ir — /il = eV to the 
system by applying /ir = eE/2 and fiR = — eE/2. We 
also rewrite the width of each miniband, = P-\-oi and 
= P — a where a, with \a\ < (3^ determines the asym¬ 
metry between the left and right energy width, so the 
relative thickness of the left and right miniband width is 
determined by a. 

The chemical potential of the cavity jXc well as the 
temperature Tc are determined by imposing conservation 
of charge and energy, Ir E Ir = ^ and J E Jr E Jr = ^ 
where J is the heat current entering from the heat source. 
From these conservation laws, we can obtain the electric 
and heat currents through the system, as well as J and 


1. Linear response 


To linear order in the temperature difference AT and 
the bias voltage E, the net current flowing through the 
system, Ir = —Ir = /, is given by 


I = GcVE GcStAT. ( 8 ) 


The electrical conductance Ge and thermopower (or See- 
beck coefficient) St of the system are 


Ge = - 


e‘^U2DA GriGri 


27rh Gri + Gri 


Q _ kB 
— — 


Gr2 E Gr3 Gr2 E Gr3 


Gri 


Gri 


(9) 

(10) 
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FIG. 2. (a), (b) Maximum power in units of 

1 ^ 20 MkET) ^ linear response with respect to 

the centers of each miniband for a symmetric configuration, 
a = 0, when (a) is /3 = ksT and (b) is /3 = bhuT. (c), (d) Ef¬ 
ficiency at maximum power normalized by Carnot efficiency 
rjc within linear response as a function of the centers of each 
miniband for symmetric setup. Panel (c)/(d) is a correspond¬ 
ing efficiency for the case of (a)/(b). All plots are obtained 
for T = 300iF and AT = IK. 


where we have introduced the auxiliary functions 


G 


rl 


Gr2 


Gr3 



dE, 


1 



dE, 


EjkBT 

l + e^GkBT' 





( 11 ) 

( 12 ) 

(13) 


The electrical conductance shows that Gli{Gri) is pro¬ 
portional to the electrical conductance of the left (right) 
superlattice, so the net conductance Ge is simply the 
series combination of the two conductors. The three- 
terminal thermopower St is determined by the differ¬ 
ence between left and right two-terminal thermopower of 
each superlattice, and vanishes when they are identical 
[laiiH]. Therefore, depending on the magnitude of the 
left and right properties, St shows whether the system is 
analogous to a p-type or n-type semiconductor. When 
St is positive for {Gl 2 + Gl3)/Gli > {GR 2 + Gm)/GRi, 
the system acts as if the mobile charge carrier is positive 
and behaves like a p-type semiconductor, and vice versa. 

The bias voltage V applied against heat driven charge 
current generates a finite output power P = \IV\. The 
power vanishes either when no bias voltage is applied or 
at the stopping voltage, Vstop = —5'tAT. The output 
power takes its maximum value at half of the stopping 


voltage. 


Pma. = \Ge\^-^^^. (14) 

In order to evaluate the efficiency 77 given by the ratio of 
output power to input heat current, we need to find the 
heat current J = —Jl — Jr injected from the heat bath. 


J = GJiVp {GeStli + Hti + Ht2)AT, (15) 


which shows the Peltier effect. Here the coefficients H, 
and Ht 2 are defined as 


n = 

Hti = 


Gr2 + Grs Gl2 + Gl3 


Gri 


2 — Gri 


Ht2 = 


ksT 
e 

27rfi ^ Gri + Gri 
{Gl2 + Gr^Y 
Gri 


Gri 
{Gr2 + Grs^ 
Gri 


(16) 


^Gr2 + Gr3){Gr2 + Grs) 


klT[{GR^P2GR,-GRe) 
-\-{Gr4 + 2Grs — Grq)] , 


together with the auxiliary functions 

{E^lkRTf 

\ _j_ ^EzIksT ’ 

E 


GrA — / JEz 
JEl- 

Grb = / dEz 
JEl- 


kuT 




nEl+ 

Gr6 = / dE.Lhi-e-^^/'^^^). 

Je:- 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


Rewriting the heat current in Eq. (15) as J = HI 


{Hti + i 7 t 2 )AT, we see the meaning of the coefficients 
more clearly. The first part is simply proportional to 
the charge current, and H shows the presence of heating 
from the Peltier contribution. The relation H = —TSt 
from Eq. ([To|) and Eq. (16) shows the coefficient H can 


be considered as the back action counterpart to St, and 
the minus sign of the relation comes because J is the 
heat current injected from the heat bath into the cavity. 
This relation shows our system satisfies Onsager sym¬ 
metry resulting from the time reversibility, which relates 
the Seeback and Peltier coefficients [49] . The second part 
shows the heat generated by the temperature difference 
and the thermal conductance is Hti +i 7 t 2 - The first one, 
Hti, contributes to the heat flow from the asymmetric 
superlattices and disappears when the left and right su¬ 
perlattices are symmetric. 

The heat current at maximum power takes the form 


J. 


maxP — 


GeStU 


Hti + H 


t2 


AT. 


( 22 ) 
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FIG. 3. (a), (b) Maximum power in units of 

1 ^ 20 •MksT) ^ kB^T y linear response as a function of 

the center of right miniband 8^ and the asymmetric param¬ 
eter a/yd, when (a) is [3 = ksT and Sq = —AksT, and (b) 
is /3 = bksT and Sq = —20kBT. (c), (d) Efficiency at maxi¬ 
mum power normalized by Carnot efficiency r]c within linear 
response as a function of 8^ and a/^. Panel (c)/(d) is a 
corresponding efficiency for the case of (a) / (b). All plots are 
obtained for T = 300iF and AT = IK. 
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Therefore, the efficiency at maximum power is given by 


'HmaxP 


GeSf 


2GeStU + AHti + AHt2 
GeS^T 


2GeStU + AHti + AHt2 


AT 


Vc, 


(23) 


where the second approximation comes from the Carnot 
efficiency for the small temperature difference, rjc = 
AT/Tc ~ AT/T. Therefore, the combination of the coef¬ 
ficients, inside of the vertical bars in Eq. (23) determines 
the efficiency ratio to Carnot efficiency. 

Fig. shows the maximum output power and the cor¬ 
responding efficiency on equal size bands {a = 0 ) as a 
function of the centers of minibands of the two superlat¬ 
tices, 8q^^ . Both the power and efficiency are symmetric 
with respect to an exchange of 8q and 8^. On the one 
hand, the maximum power arises when one of the two 
center value stays around twice of p and the other center 
is deep below the equilibrium chemical potential, below 
about —3y5. The center around 2/3 is the position that 
makes the bottom of the miniband stay around the equi¬ 
librium chemical potential, and the other center which 
is below —3/3 makes the top of the miniband deep be¬ 
low the chemical potential, because the miniband width 
is 4/3. On the other hand, the efficiency acts symmetri¬ 
cally depending on the position of band centers as well 
as the power. However, maximum efficiency comes when 
one of the center of miniband is bigger than 2/3 with an 



FIG. 4. The figure shows nonlinear response when the mini¬ 
band width \s P — ksT. (a) the maximum output power as 
a function of temperature difference AT/T for the optimized 
parameters 8o = —^ksT, 8 q = 2/cbT, and a/ /3 — —0.25. (b) 
shows the corresponding efficiency. The power is in the units 

^2D^{hBT)^ 

2TTh 


appropriate value of the other center, and it comes with 
suppressed output power. Therefore, depending on what 
we want to optimize, the output power or efficiency, we 
can chose the position of band centers. 

We also show the maximum power for the different 
miniband widths from the top panels of Fig. As the 
electrons transport only within the miniband, a wider 
miniband allows more electrons to transport and generate 
more power. However, as the miniband width goes too far 
above ksT, the power increase stops. The reason is that 
the energy window of the order ksTf^jj. will be a more 
effective energy guard for the carriers. Different from 
the power, efficiency is reduced as the miniband width 
increases, Eq. (i c),(d)). As the width increases, the 
energy filtering by the superlattices is lesser efficient so 
the efficiency decreases. However, as the width increases 
continuously, the efficiency will saturate at some point 
with the same reason for the power saturation. 

These results show that the superlattice heat engine 
has more power output than the quantum well engine 
m- The maximal output power of the superlattice 
heat engine with units Po = iy 2 DAkBT/{2h{kB^T/2)‘^) is 
Pmax ~ 1 - 8^0 with the efficiency r]maxP - 0.05r]c (for the 
miniband width, /3 = kBT), while the quantum well heat 
engine is Pmax ~ with the efficiency r]maxP - O.OTr^c* 
From Eq. (nnH’ we can give a simple reason for this: 
when the miniband width is suppressed (longitudinal en¬ 
ergy window becomes one value and the electrons trans¬ 
port only with the certain longitudinal energy), our re¬ 
sults approach the quantum well case. Therefore, a large 
miniband width permits more electrons to transport in 
the longitudinal degree of freedom. While the superlat¬ 
tice engine is more powerful, it is less efficient an energy 
filter than the quantum well case: this is from the differ¬ 
ence between the miniband of the superlattice and the 
sub-band threshold of the quantum well. 

We now turn to an asymmetric system, a ^ 0. To 
find the optimized combination of the centers of mini¬ 
band and the asymmetric parameter a for the maximum 
output power, we consider the output power and the effi- 
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ciency ioi Fig. shows plots of the power as a 

function of a and £q. When /3 = ksT^ the largest output 
power Pmax ~ 1.85Po (around a 2% increases from the 
symmetric case) arises for a ~ —0.2/3 with £^ « 2 . 6 /cbT. 
In this case, the efficiency goes to r]maxP ~ O.OGr^c (a 20% 
increase from the symmetric case). Meanwhile, the max¬ 
imum efficiency 77 ~ 0.16?7c comes with a —/3 and 
£^ bksT for P = ksT. When a the miniband 

widths become /3^ —> 0 and Pr 2^5, and the output 
power is strongly suppressed. We can understand the 
reason for this last fact because when the width of the 
left miniband vanishes, transport from the left is cut off, 
giving no power, but very good energy filtering, which in¬ 
creases efficiency. To explain other features of the plots, 
we note that for the right superlattice, the center of mini¬ 
band stays slightly above the equilibrium chemical po¬ 
tential, and the width of miniband A = A^r is larger 
than the width of the distribution function ^ ksT. Con¬ 
sequently, making the miniband wider has no effect on 
the electron transport because the number of electrons 
with energies around minimum/maximum of the energy 
window is exponentially small. Simultaneously, for the 
smaller miniband of the left barrier, centering it around 
the energy region where the left reservoir is occupied but 
the cavity not, gives very good current production. As 
the miniband width [3 increases. Fig. &) , a moves to 
a positive value while the center £^ stays around the 
value that maximizes the power for the symmetric case. 
Therefore, depending on the size of the miniband width 
relative to the width of the occupation function of the 
reservoir and the cavity, a can be fixed to give the opti¬ 
mized results. 


2. Nonlinear response 


It is interesting next to consider the output power and 
the efficiency in the nonlinear regime, where qualitatively 
new physics can appear [5Qll52] . We numerically calcu¬ 
lated the stopping voltage Vstop^ the centers of miniband 

£q^^ ^ and asymmetry parameter a in order to maximize 
the output power. For the case p = ksT^ we show the 
maximum output power and the efficiency in Fig.[^ This 
is obtained for the miniband centers £q = —^ksT and 
£^ = 2k rT with the asymmetric parameter a = —0.25/3. 
The output power increases quadratically in the temper¬ 
ature difference for the fixed average temperature T as 
in the linear case, Eq. (14), but it also depends on the 
average temperature, Pmax r\j T AT^. Therefore higher 
average temperature as well as the temperature difference 
gives bigger output power. 

In Table |lj we compare the maximum output power 
and the corresponding efficiency for three systems: the 
quantum dot [9], quantum well m, and superlattice 
based three terminal heat engines for a realistic device 
parameter rrieff = 0.067me with the room temperature 
T = SOOiC and the temperature difference AT = IK. We 
see that Pmax of the superlattice heat engine generates 



Quantum Dots 

Quantum Wells 

Superlattices 

Pmaxi'W/cm'^) 

0.1 

0.18 

0.3 

PmaxP {j]c) 

0.2 

0.07 

0.06 


TABLE I. We compare our result with other three terminal 
geometries: (power-optimized resonant width) quantum dot 
and quantum well heat engines for T = 300A and AT — IK. 


a larger power about 0.3W/cm^ with a small reduction 
for the efficiency O.OOr^c* Therefore, the superlattice heat 
engine is the more powerful heat engine. 


C. Quantum refrigerator based on superlattices 


Now, we consider the same geometry but for a differ¬ 
ent purpose, a cooling system. If proper bias voltage eV 
is applied over the junction and the minibands of super¬ 
lattices are suitably arranged, a current flows from right 
to left as hot electrons tunnel through the miniband of 
the left superlattice from the cavity to left reservoir and 
cold electrons from the right reservoir tunnel through the 
right miniband to the cavity. This leads to decrease in 
the average energy of electrons in the central cavity, that 
is cooling utilizing the Peltier effect: in the Peltier effect, 
the junction is electrically biased and a produced heat 
current flow given an electric current. 

Eor the refrigerator, AT is negative in our notation, 
and the temperature of the center cavity cooled down to 
the amount of |AT|. The purpose of refrigeration is to 
achieve a large temperature reduction. The base temper¬ 
ature of the refrigerator is defined as the temperature for 
equilibrium where the evacuated heat current balances 
any external heat leaks. Since we are now considering no 
external heat leaks, such as electron-phonon coupling, 
the base temperature is the temperature for which J = 0 
in Eq. (15) with the temperature reduction 


ATo = - 


an 


GeStU + Hn + H, 


-E, 


T-t2 


(24) 


in the linear regime. The refrigeration works only when 
the heat current is emitted by the cavity into the reser¬ 
voirs, J > 0, therefore the temperature reduction of the 
cavity should be in the range of ATq < AT < 0. 

The applied bias voltage V and the absorbed heat from 
the cold cavity let the heat flow J rejected to the left and 
right reservoirs. Therefore, the cooling power is J, and 
an efficiency of the cooling is normally characterized by 
the coefficient of performance (COP), defined as the ratio 
of the cooling power to the total input power P = IE, 


(/) = 


J 

P‘ 


(25) 


Similar to the efficiency of the energy harvester, COP is 
also bounded by the Carnot value, 0 < = (pc- 

Moving beyond the linear regime is important to find 
the nonlinear behavior of the refrigerator, which deter¬ 
mines the lowest temperature it can reach. Eig. Ha) 
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AT/T 


(b) 



0.012 
0.006 ^ 
0 


FIG. 5. (a) The nonlinear cooling power is shown as a function 
of the temperature difference and the bias voltage when the 
other parameters, Sq , 8^, yd, and a, are optimized to give the 
maximum cooling power. The cooling power is in the units of 
i^ 2 D^{kBT) ^ shows the heat current and the COP 

as a function of the temperature difference when the bias is 
chosen to give the maximum cooling power, eV ~ 2k bT. 


shows a plot of the cooling power vs applied bias and tem¬ 
perature differences, when other variables, 8^ = — 

8q = 2 /cbT, [3 = ksT and a = 0.3/3, are optimized to 
give the largest cooling power for the positive applied 
bias y > 0. Here, a > 0 means that left miniband width 
is wider than the right miniband, and the reason can be 
understood as the energy harvesting that we discussed 
Note that, the position of the left/right 


IIB 


in section _ 

miniband of the refrigerator is opposite from the energy 
harvesting, so this makes the conduction electrons take 
energy from the cavity and continue to cool it. 

The black curve in panel (a) of Fig. (§ represents when 
the cooling power becomes zero, and the right hand side 
of the line (within the parabola) is for the positive cool¬ 
ing power which is the region that works for the refrig¬ 
eration. Therefore, we can see the optimal region for 
the temperature difference and the applied voltage, as 


shown in the plot. When the temperature difference is 
anywhere between the zero and ATmax^ as an example 
the vertical line in the figure, the applied voltage should 
stay in between Vmin and Vmax to make the engine work 
as the refrigerator. Specially, when the applied voltage 
is around 2k bT^ the refrigerator gives maximum cooling 
power for a given temperature difference, and we also 
can have the maxumum temperature reduction |ATo| of 
the system. When |AT| 0, the heat current increases 
linearly in panel (b). So, the maximum temperature re¬ 
duction for this bias voltage is |ATo| = \ATmax\ — 0.1 IT. 
For example, at room temperature, the maximum tem¬ 
perature reduction is \ATmax\ = 30iF by applying the 
bias around 50meV. If the temperature reduction is 
|ATo| = 0.05T = 15iF, the corresponding cooling power 
is approximately llkW/cm^. The COP at eV = 2 /cbT 
as a function of the temperature difference is also plotted 
in Fig. ([^b)). The maximum COP <p ^ 0.015(/)c comes 
with when | AT| ^ 0.06T. Therefore, the refrigerator can 
be operated at the optimal regime for the cooling power 
or COP by reasonably choosing the parameters. 


III. HEAT TRANSPORT BY PHONONS 
A. Reduced phonon heat current 

One important consideration in thermoelectric heat 
engine is the heat fiow carried not only by the conduction 
electrons, but by the phonons as well. Here, we calcu¬ 
late this effect for our system to see how it will affect the 
efficiency in a more realistic modeling. 

The major heat fiow from hot to cold reservoirs is car¬ 
ried by excitations such as phonons. For this three ter¬ 
minal heat engine, the phonon heat fiow is in parallel 
with electronic fiow J®, and the total heat fiow for a given 
generated power P is the sum of them, J = 

Then we can rewrite the efficiency in the present of the 
phonons 

The phonon heat current in the superlattice differs 
from the heat current based on the bulk material prop¬ 
erties due to the new periodicity of the structure. The 
presence of interfaces can alter the phonon spectra re¬ 
sult from wave interference scattered at the interfaces. 
The formation of miniband gaps in superlattices due to 
phonon interference leads to a reduction of the phonon 
velocity Vg{k) which gives a reduction of the thermal 
conductivity. To illustrate this effect in a simple model, 
we assume the complete separation of longitudinal and 
transverse vibrations for phonons with the wave vector 
to parallel to ^ axis. Also assuming perfect interfaces, 
the transverse momentum qx{y) is conserved. When we 
consider phonons propagating in the growth direction, 
this can be treated from a Kronig-Penney type approach 
for one dimensional atomic chain [34]. We suppose the 

















(a) OJ 


(b) To 



FIG. 6. The phonon dispersion curves from a one dimensional 
atomic chain for the case of mass ratio (5 = 2.6 in the extended 
zone representation, (a) and (b) show ni = n 2 = 1 and 
ni = 712 = 3, respectively. The dimensionless parameters are 
dehned as q = qa and UJ = uj^/Mi/g. 


same monolayer spacings a and the magnitude of atomic 
constants g in between all the atoms. The layer one (two) 
has ni{n 2 ) atoms with mass Mi (M 2 = ^Mi), and the 
thickness of sublattice one (two) is di = riia (^2 = n 2 a) 
which gives the length of unit period d = + (i 2 . Then 
the characteristic equation is 


cos{qzd) = cos{kidi) cos 
1 — cos{kia) cos(/c 2 a) 
sm{kia) sin(/c2a) 


(fe<^2) 

sm{kidi) sin(/c2C^2), 


(27) 


where /ci( 2 ) is the phonon wave vector of each layer 1 ( 2 ) 
with cos(/ci( 2 )a) = 1 — Mi( 2 )<^^/ 25 f, and qz is a longitu¬ 
dinal superlattice wave vector. This model gives us the 
dispersion of the longitudinally polarized phonons for the 
cross plane transport with zero transverse momentum. 
The dispersions for different periods are shown in Fig. 
for 771 = 772 = (I 5 3) with the mass ratio S = 2.6. Increas¬ 
ing the superlattice periods (771 and 772 ) gives more band 
folding and decreases the average group velocity. 

At nonzero transverse momentum, the longitudinally 
polarized phonon mode for the cross plane dynamics is 
described as [36] 


^(2 - cos(gj,a) - cos(g^a)), (28) 

where ujq^ is a solution from Eq. ^7\ and fit is a charac¬ 
teristic frequency of the material. Therefore, the disper¬ 
sion of non zero transverse momentum can be used to cal¬ 
culate the group velocity in the super lattices, Vg = V^cjq. 
Since we are dealing with the heat transfer only to the 
cross plane direction, the group velocity that we need is 
only the ^ direction, 

Vg,z = —Vg,ld, (29) 

CJq 

where we use the one dimensional group velocity Vg^id = 
dujg^jdqz. Now, we can write the phonon heat current 
from the reservoir r into the cavity through the superlat¬ 
tice, 

= I ^ hw^Vg^MTiqz) ABiq,Tr,T,),im 


where AB{q,Tr,Tc) = br{^ — bc{^ is the occu¬ 
pation difference between the reservoir and cavity 
with the Bose-Einstein distribution function brjdg) = 
{ex.p{bw{q)/kBTr/c) — 1)~^ and T{qz) represents the 
transmission function which depends only on the longi¬ 
tudinal momentum qz. In our analysis, we assume the 
ballistic case with T{qz) — 1- When we put the longi¬ 
tudinal group velocity Eq.(|2^ into the first equation of 
Eq. ( |M| ), we can factorize the heat current to the longi¬ 
tudinal (z) and transverse (x, y) directions 

VgM<lz) A7’(9., Tr, T,), (31) 

with a modified occupation difference 


AF{q^,Tr,Tc 



dqxdqy 

{2ny 


AB{q,Tr,Tc). (32) 


This function gives the effective phonon occupation dif¬ 
ference per unit area between the cavity and reservoir due 
to the transverse momentum. The transverse momen¬ 
tum of phonons contributes only through the distribution 
function, and the energy and velocity dependencies of the 
heat current are only from the longitudinal momentum 
qz. This analysis suggests that the heat current can be 
treated as a one dimensional calculation with an effective 
occupation function. We rewrite the new expression in 
terms of an integral over longitudinal frequency instead 
of the momentum space 


ni+n2 


= T. J ' AF{lU,^ , Tr, B) 


i=l 
ni+n2 

^ /■ 

co¬ 


ded, 


Qz 


277 


hujg^AF{ujg^,Tr,Tc), 


(33) 


where uj- is the minimum (maximum) frequency of 
the 7^^ miniband of longitudinal momentum from Eq. 
^7\ . The conversion to frequency space introduces the 
phonon density of state of the superlattice, and we have 
= {vg,id7T)~^• Therefore, the group velocity 
and the density of states cancel out and we have a simple 
equation of the phonon heat current in Eq. (33). 

When we consider for the electron transport of the su¬ 
perlattice, note that only the lowest miniband was con¬ 
sidered, because the energy gap of the first and second 
miniband is larger than the thermal energy ksT^ for the 
typical materials of the superlattice. So, when the chem¬ 
ical potentials of the reservoir and the cavity stay around 
the first miniband, energy window ksT of the occupation 
will guide the transport only through the first miniband. 
Therefore, it is a reasonable approximation to ignore the 
higher minibands. In contrast, the Bose-Einstein distri¬ 
bution function in the frequency domain is broader than 
the maximum frequency of the acoustic dispersion, and so 
all of the acoustic phonons contribute to the heat trans¬ 
port. 
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(a) 

/P*'xlO'‘[W/cm2] 



(b) 

JxlO'‘[W/cm^] 



FIG. 7. The plots show electron and phonon heat currents of 
GaAs/AlAs super lattice for the case of the maximum output 
power with optimized parameters with AT = IK. (a) shows 
the heat current by phonons when the atomic layers are 
ni = 712 — (I 52 ), with the bulk AlAs case, (b) shows the 
magnitude comparison between K and . 


The total phonon heat current from the heat engine 
is also obtained from the heat conservation + 

= 0. However, different from the heat current by 
the electrons, the difference of the phonon heat current 
through the left/right reservoir is only determined by the 
temperature difference between the cavity and reservoirs 
which is contained only in the modified occupation differ¬ 
ence AF{ujq^,Tr,Tc). The temperatures of the left/right 
reservoir are the same in our case, therefore, the phonon 
heat current conservation leads to the total phonon heat 
current such as only if there is symmetric 

heat conductance. 

In order to calculate the effect of the phonons on the 
thermodynamic efficiency of this engine, we consider a 
low temperature regime, where an analytic investigation 
can be made, followed by a numerical investigation of the 
high temperature regime. 

In the low temperature case, we assume the tem¬ 
perature only excites low energy modes. Up to a 
small numerical factor, this permits us to approxi¬ 
mate the Bose-Einstein as a Boltzmann distribution, as 


(a) rf/Tjc 



(b) nhc 



FIG. 8. Efficiency at maximum power of GaAs/AlAs super¬ 
lattices for the value of AT = liF, /3 = ksT, and ni = 77-2 = 2, 
and are chosen to maximize power. The efficiencies by 

the electron heat current (a) and the total heat current (b) 
are shown as a function of average temperature. 


well as only consider the long wavelength modes, so 
the dispersion relation may be approximated as, cj ~ 

'\J^qz + ffi Q-y)/^' Computing the effective dis¬ 

tribution by integrating over Qx^Qy^ we find. 


F{uJqz) 


hWqzkBT hojgJkBT 

27ra^(fiUt/4)^ 


(34) 


This allows us to approximate the heat current per unit 
area from the first mini-band as. 


jph 


{'KaMltYh 

16 

{TTahritYh 


[{kBTLf - [ksTef] 

{keTfAT, 


(35) 

(36) 


where the last limit is in the linear response limit. This is 
consistent with a Debye treatment of the phonon trans¬ 
port. This result shows that for low temperatures, the 
phonons freeze out, and the energy is predominately car¬ 
ried by the electrons. 

Moving on to the high temperature limit, we consider 
a numerical investigation. This limit is quite impor¬ 
tant for room temperature applications in mind. Eig. 

shows the heat currents for GaAs/AlAs superlattices 
when AT = liF, as an example. The phonon heat 
currents of superlattices and bulk material in panel (a) 
increase with temperature increasing, but saturate at 
high temperature. Moreover, as the superlattice atomic 
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layer increases, the phonon heat current decreases and 
saturates when the atomic layers are ni,n 2 > 10 [38] . 
The practical superlattices usually have the atomic lay¬ 
ers ni,n 2 > 10, so, the practical phonon heat current 
of the superlattices has about 40% reduction of the heat 
current from the bulk material around the room tem¬ 
perature. Panel (b) compares the electron and phonon 
heat currents for the optimized parameters to give the 
maximum output power. As the temperature increases, 
the electron heat current increases quadratically while 
the phonon heat current saturates, so the electron heat 
current reaches the phonon heat current at high temper¬ 
ature. However, at room temperature, the phonon heat 
current is still an order of magnitude higher than the 
electron heat current. 

We compare the efficiency by the electron heat current 
and the total heat current in Fig. The phonon heat 
current dominates over the electronic one, ^ J®, 
when the temperature is around room temperature. In 
this situation, we can rewrite the total efficiency 

y^e+ph _ pjjph^ ( 37 ) 

The efficiency depends only on the phonon heat current 
and the generated power. Therefore, we have the maxi¬ 
mal efficiency when the power is maximal, and the power 
and efficiency reach maximal value together. Because the 
phonon heat current is an order of magnitude higher than 
the electron heat current, the total efficiency decreases 
about an order of magnitude. 


B. Optimized condition for minimum phonon heat 
current 

Practically, there could be a way to reduce the phonon 
heat current more than we showed based on the sim¬ 
ple theory and example. Based on the one dimensional 
atomic chain model, we have calculated the phonon heat 
current in the growth direction and estimated their con¬ 
tributions to the efficiency in the heat engine based on 
GaAs/AlAs superlattices. Our approach shows that the 
reduction of the phonon heat current comes from the re¬ 
duction of the group velocity near the folded Brillouin 
zone edges. This calculation also assumes the ballistic 
transport T{qz) — 1- However, as the constituent lay¬ 
ers become thicker than the phonon mean free path, the 
phonon transmission function needs to be modified as 
Tfe) < 1. In this case, we also need to treat the phonons 
as particles and use the theory such as the Boltzmann 
transport equation [37]. Moreover, other mechanisms 
such as phonon spectra mismatch and scattering aris¬ 
ing from the roughness of the layer interfaces also play 


an important role to understand the reduction of the ex¬ 
perimental results [38l [53ll56] . 

Some of high figure of merit thermoelectric materials 
show reduced lattice thermal conductivity in the super¬ 
lattice structures, for example Si/Ge or Bi 2 Te 3 /Sb 2 Te 3 . 
The reduction of phonon heat current can be maximized 
by the proper choice of superlattice period compared to 
the mean free path of the phonons: a theory predicts 
the thermal conductivity minimum as a function of layer 
spacing [39] , and some works show the minimum thermal 
conductivity depending on the superlattice period and a 
ratio of the layer thickness for the materials [iniiiuiiT]. 
Therefore, these materials, instead of our example, with 
optimized superlattice period will give more reduction of 
phonon heat current of the system. 

IV. CONCLUSIONS 

In this paper, we investigated a heat engine by thermo¬ 
electric effects in superlattice structures in three termi¬ 
nal geometry. First, our work considers the engine as a 
energy harvester, and shows the advantages of superlat¬ 
tice heat engine in the large output power compared to 
the other similar geometry due to the box shaped trans¬ 
mission function which comes from the electron energy 
miniband. Our theory predicts that the maximum power 
under optimized conditions can be larger than similar 
resonant tunneling devices, with comparable efficiency 
at maximum power. Second, a different regime of the 
system parameters makes the engine works as a refriger¬ 
ator, and we shows the optimized regime of the param¬ 
eters for the maximum cooling power and coefficient of 
performance. In addition, we analyzed the phonon heat 
current to find the total efficiency by the performance of 
electrons and phonons together. The reduction of phonon 
heat current in the superlattice compared to the corre¬ 
sponding bulk material offers higher total efficiency. The 
trade off between power and efficiency is overcome by 
the reduction of the phonon heat current which can be 
achieved either by operating at low temperatures, or by 
engineering the system to have low phonon conductivity, 
while keeping high electron conductivity. Easy fabrica¬ 
tion of these devices with advantages for both the power 
and efficiency show this heat engine is a promising device 
for next-generation thermoelectrics. 
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